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Abstract 

The CDF collaboration has recently reported a large deviation from the standard model of 
the ti forward-backward asymmetry in the high invariant mass region. We interpret this mea- 
surement as coming from new physics at a heavy scale A , and perform a model-independent 
analysis up to ©(l/A^) . A simple formalism to test and constrain models of new physics is 
provided. We find that a large asymmetry cannot be accommodated by heavy new physics 
that does not interfere with the standard model. We show that a smoking gun test for the 
heavy new physics hypothesis is a significant deviation from the standard model prediction 
for the ti differential cross section at large invariant mass. At M^i > 1 TeV the cross section 
is predicted to be at least twice that of the SM at the Tevatron, and for M^f > 1.5 TeV at 
least three times larger than the SM at the LHC. 

1 Introduction 

The top quark is the most massive point-like particle known to exist. As a consequence, within 
the Standard Model (SM), the top is largely responsible for the hierarchy problem. Furthermore, 
in most natural models it is linked to electroweak symmetry breaking. Therefore, there is strong 
motivation to search for new physics effects associated with top physics. 

The CDF collaboration has recently announced several intriguing new measurements that ex- 
hibit large deviations from the corresponding SM predictions. Evidence for an anomalous forward- 
backward tt production asymmetry was observed for large invariant mass of the ti system [1]: 

A%f^ = A"(Mit- > 450 GeV) = +0.475 ± 0.114 , (1) 

to be compared with the SM prediction [21 El g], A%q = +0.09 ± 0.01. Previous DO and CDF 
measurements of the inclusive ti asymmetry [5l [6] also show deviation from the SM prediction. 
Another recent CDF analysis in the dilepton channel [7] supports this deviation, and furthermore 
finds a rising Mu dependence for the forward-backward asymmetry. 

Additionally, the CDF collaboration has recently made progress in studying the mass distri- 
bution of highly boosted jets (pr > 400 GeV for the leading jet) |8], and found a hint for an excess 
of events in the high mass region [9]. 

The above measurements suggest that new physics affecting the top sector is present. Our 
approach in this work is the following. We interpret the measurement of A45Q in terms of new 
physics, checking the consistency of such a scenario with other measurements that do not show 
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any significant deviation from tlie SM predictions. We furtlier discuss tlie effects of sucli new 
pfiysics on ultra-massive boosted jets at tlie Tevatron. We then make predictions for the invariant 
mass distributions of top pairs soon to be measured at the Tevatron and LHC. 

Several works have interpreted the recent CDF measurement of ^4450 within specific models 
of new physics [IUl[IIllI21[I31[Il[I3[inilIZl[IHl[ISll2ni[nil221[23j. Similarly, model-independent 
analyses were performed [211 ESI [20] and new physics models were invoked [271 ESI |2Hl EDI ED ESI 
EalEllEHlEniETlEElEHlllQlIlT] to explain earlier DO and CDF measurements of the inclusive 
asymmetry [3 [6] . 

We focus on the class of models in which the scale of new physics is well above the scale Mti 
relevant to the CDF measurements. The effects of such new physics can then be described from 
a low energy model-independent perspective, using the language of effective field theory. Ref. [12] 
performed a similar analysis, further assuming that the dominant contribution to the forward- 
backward asymmetry comes from interference between the new physics and SM contributions to 
top pair production. Denoting the scale of new physics by A, Ref. [12] found that in the presence of 
an axial octet operator producing a pair of tops from a pair of up quarks at 0(1/A^) , the observed 
tt forward-backward asymmetry can be accounted for. Here we relax the assumption of interference 
and consider all operators contributing to ti production up to order 1/A^. We provide a simple 
formalism that enables one to easily obtain constraints and predictions for models consistent with 
our framework. We derive model-independent predictions regarding near future measurements 
that will sharply test our general underlying assumptions. 

The paper is organized as follows. In Section |2] we review the data relevant to our study. 
Section [3] defines the set of operators in our effective Lagrangian. Section |4] relates the operators 
to the observables. Our results are presented in Section [5j In Section |6] we discuss predictions for 
hard top physics at the Tevatron and LHC. We conclude in Section [7} 

2 Relevant Data 

In this work we analyze the effect of heavy new physics on the forward-backward asymmetry at 
large Mti. Roughly, we aim to account for a new physics contribution of 

= +0.40 ±0.11, (2) 

assuming that the rest of the asymmetry in Eq. ([T]) comes from the SM. 

Other top-related measurements do not show significant deviations from the SM predictions. 
Consequently, they provide constraints on the parameter space of the effective Lagrangian. The 
first such observable is the ti differential cross section, which we choose to represent by the following 
large Mti bin [33] 

cr7oo = (t""(700 GeV < Mu < 800 GeV) = 80 ± 37 fb , (3) 

as in [12]. This is consistent with the SM prediction [2 El], cttoo = 80 ± 8 fb. In [12], the 
inclusive ti cross section was also used as a constraint. However, the theoretical estimation of 
the cross section originating from threshold effects is still under investigation (compare [331 HHl IS] 
with [11]). Furthermore, the dynamics of our heavy new physics naturally affects the measurement 
at large invariant masses more significantly. Thus in our study we do not use the inclusive ti cross 
section to constrain our parameter space. (Note though that our results below are within the 
combined theoretical and experimental uncertainties for the inclusive observables.) The same 
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argument leads us to refrain from considering A" in the low invariant mass region, as well as the 
inclusive asymmetry. We therefore use a theoretically-cleaner observable, also relevant for A^jq , 
related to the cross section above 450 GeV, (T450 . Note that there seems to be some discrepancy 
in |43j between the measurement in this range and the SM prediction. However, in a more recent 
measurement reported in ^ (but not translated to the partonic level) this discrepancy is not 
present. We thus assume that the SM prediction agrees with the measured value of the cross 
section above 450 GeV. For concreteness, we use the relative uncertainty reported in [T] for the 
450-500 GeV Mti bin (see below), which dominates the uncertainty in the Mti > 450 GeV range. 

In order to minimize the impact of next to leading order (NLO) corrections to the new physics 
(NP) contributions, we normalize the latter to the SM one. We assume that the ii'-factors are 
universal, so that the NP/SM ratios at leading order and next to leading order are the same. This 
assumption is reasonable for the effective operators generated in the SM after the highly virtual 
gluon is integrated out Additional NLO corrections are formally down by Olas/n) and are 
henceforth neglected. In practice, log-enhancements of these contributions might be present which 
could modify the analysis. A full NLO computation of the effective theory is yet to be done, and 
is beyond the scope of this work. However, as a self consistency check of our analysis, we find 
below that the SM-like operators indeed account for the dominant part of the forward-backward 
asymmetry (see Fig. |2|, supporting the above assumption. 

Combining in quadrature the experimental and theoretical uncertainties, we represent Eq. ^ 
and the uncertainty on the cross section above 450 GeV as follows: 

lA^rool = Ko/^^^ol < 0.5 , |A^45o| = KV^!^o\ < 0-2 . (4) 

It is also intriguing to explore the implications of the new physics in the context of the CDF 
boosted jets study [Bl E]. The cross section for ultra-massive boosted jets (not coming from QCD 
events) can be estimated as follows |12] 

[21-(8.7±3.1)i?^Ls] fb, (5) 

where ab is the cross section of hadronically-decaying ti with a pt cut of 400 GeV on the leading 
jet and Rmass is a parameter that determines the QCD background, as defined in |l2l HE]. (An 
assumption of naive factorization of the jet mass distribution yields -Rmass = 1 , while matched 
Monte-Carlo simulations give -Rmass = 0.87 [l2l HH] with an excellent agreement on this value 
between the different generators.) The SM prediction for the top contribution is af^ = 2.0 ± 
0.2 fb [19]. We interpret the excess as top pairs, generated by the new physics source. The 
magnitude of this effect is then ^2] 

m ^ ^r/crr = 5 ± 2 , (6) 
where af^ is the new physics contribution to the boosted cross section, assuming -Rmass = 0.87. 

3 The Operator Basis 

As stated above, the basic assumption that we employ is that the new physics is characterized by 
a scale A that is larger than the invariant mass of the top pair M^f in the measurements which 
we consider. The natural approach is then to use a set of effective operators to describe the new 
physics. These operators must lead from an initial uu state to a final tt state, and as such appear 
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at dimension six and higher. (The contribution of dd — > ti at the Tevatron is at most 15% that of 
uu — 7- ti for Mft above 450 GeV, as relevant for the observables that we consider.) 
At (9(1/A^), there are only two four-quark operators that interfere with the SM: 

o^y = (u^^n) (trT^t) , o\ = {uj^^'n) (tri'T-t) , (7) 

where the superscript 8 denotes an octet color structure. Allowing for contributions that do not 
interfere with the SM, there are two more vector octet operators at this order: 

0\v = {u7,7'T^u) {t^^TH) , OIj, = {u^,n) {tn'T'^t) . (8) 

There are four additional orthogonal combinations of color contractions, given by: 

The list of dimension six operators is concluded with eight scalar and two tensor operators: 

Ol'^ = (uTi^su) (tTi,8t) , O'/ = (uTi^sl'u) {tT.^sl't) , 

= I [u T^^su) {tT,^8l't) , 0];l = t {u T.^sl'u) {m,st) , (10) 
O'/ = {uTi,sCT^''u) (tTi,8a^,t) , 

with Ti = 1 and Tg = T'^ . 

The above dimension six operators contribute to top pair production at (9(1/A^) as well, via 
the square of their amplitudes. Another type of contribution at (9(1/A^) comes from chirality- 
conserving dimension eight operators that interfere with the SM. These can be constructed by 
applying two covariant derivatives in various ways to the operators in Eq. ([T]). However, naive 
dimensional analysis shows that their value is given by c^/(167r^) , where c is the coefficient of any 
dimension six operator. This condition generically holds in case of strongly-coupled new physics, 
i.e. when the NP scale is roughly 5-10 TeV |12]. Moreover, even in models with a lower scale {e.g. 
a ~2 TeV s-channel resonance), producing a large value for ^4450 typically leads to a suppression 
of the dimension eight contributions compared to the square of the dimension six amplitudes. We 
thus neglect these dimension eight contributions in what follows. 

Note that in principle there are also dimension six chromo-magnetic/electric u and t dipole 
operators that can be considered. Their effects at 0(1/A^) in the hard regime were shown to 
be negligible in |12]- As they involve chirality flips, their contributions at order 1/A^ are suppressed 
by at least (mt/A) compared to their 1/A^ effects. There are also chirality-flipping dimension eight 
operators which interfere with the SM, and are again suppressed by the same factor. 

To conclude, we describe the hard region of the tt physics by the following effective Lagrangian: 

i 

where the real coefficients and the operators Oi are hsted in Eqs. (|7|)-(|l0|)- Below for 

simplicity of notation, q will denote Cj/A^gy, where Axev = A/TeV. 

In our analysis we perform all calculations at leading order and neglect renormalization group 
running and mixing. Consequently, we also do not discuss the contribution from operator mixing 
to dijet production at the LHC [13]. In principle, the NP can couple to light and heavy quarks 
with different strengths. Generically, the operators discussed above can be characterized by a 
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mixed coupling of NP to light and heavy quarks, gun x gu, whereas the constraints from dijets are 
sensitive to operators characterized by a strength gf^^. Examining the present bounds from dijet 
production at the LHC [501 EIIj a hierarchy of guu/gtt ~ 1/5 is required in order to comply with 
the data. 



4 Relating Operators to Data 

We now write the contribution of the operators to the various observables of interest. We first 
focus on the vector operators of Eqs. ([Tjl-Q. 

It is natural within the vector sector to distinguish between the operators that interfere with 
the SM and those that do not. The latter set of operators can be parameterized by: 



^ 2 



(12) 



I { {4a ± c\yY + I [(4 ±C\Y+ ± c\yY] } , 

= w'^ + , tanO = W-/w-f- . 
The relevant observables of Sec. [2] then take the simple form 

Nx ^ ax4 + bxicl.y + dxic\y + exR\ (13) 
Ago = («4 + + cos 29^ {1 + iV45o)"' , 



(14) 



where the subscript X = 450, 750, b and in Nx we neglect a term which is proportional to sin 26 
and suppressed by 4mj/M^|. The coefficients {a,b,d,e)x and a, (3 are pure kinematical factors 
given by 

(a, 6, rf, e)45o = 0.35,0.043,0.023,0.033, (15) 
{a,b,d, e)7oo = 0.76,0.16,0.11,0.14, (16) 
{a,b,d,e)b = 1.5,0.57,0.46,0.51, (17) 
a,(3 = 0.17,0.043, (18) 

where we use the MSTW parton distribution functions [52j at leading order in this calculation. 
The physical interpretation of R is very clear — it parameterizes the overall size of the operators 
which do not interfere with the SM. The angle 6 controls how much these operators project on the 
asymmetry. For a given R, the asymmetry is maximized for cos 2^ = 1, justifying the omission of 
the sin 26 term above. 

It is useful to obtain relations between the various observables, allowing for a simple estimation 
of the new physics contributions in terms of the constraints. One such relation is between the 
boosted tops enhancement factor and the two cross section constraints, given by 

iVb = -0.12 4 + 5.6 iV7oo- 7.5 iV45o, (19) 

such that 



4 + 10 < 10Vl + iV450 - 0.21 iVroo and 4 + 21 > 15.6 V1.8 + iV450 - 0.24 A/'too . (20) 



The constraints in Eq. (20) define the range in which R>0 and is real. 
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Another useful relation is between A^jq and the constraints A^45o and A^roo- To obtain this, we 
can substitute c\ for A^7oo using Eq. ([T3|, yielding 



{I + N450) X A%Q= (0.51 + 0.13 4) ^JN^oo- [0.76 4 + 0.16 (4)2 + 0.14 + 0.022 cos 2^ . 

(21) 

This relation is only valid for c\ > , and indeed, as we show below, accounting for A^jq as in 
Eq. (|2| requires c\ > . 

We would now like to comment on the most general case where the scalar and tensor operators 
of Eq. (10) are also present. It is straightforward to show that the effects of the latter can be 
captured by our Eqs. (13)-(18) with the following redefinitions 



(22) 
(23) 



COS26' F(6',V^,6'y,6'5T) = cos26'cos^6'y + y-cos2^/'sin^6'ycos^6'5T. 
i?2 ^ujI + wI+ rlj, + rp , 

with |F| < 1, and where we defined 



ST 



2 I 2 

y+ + y- 



+ 



y± = yc^T ± 



and 



4± 



'3 1 



(24) 
(25) 



+ ^2 = R COS 6v , r ST = Rs^T^(^v COS 6 ST , rp = RsinOv sinOsT , tan?/^ = . (26) 

The physical interpretation of the additional parameters is as follows. rsT,p represent the overall 
size of the scalar/tensor and pseudo-scalar operators; R being now the overall size of all of the 
non-interfering operators. The angles 6v,st determine the distribution of R among the vector, 
scalar/tensor and pseudo-scalar operators, while the angle 6 (ip) parameterizes how much the 
vector (scalar/tensor) operators project on the asymmetry. Note that in the case where vector 
operators are absent, one has |F| < v/273. 



5 The Forward- Backward Asymmetry 

The stage is now set to study the parameter space that explains the large forward-backward ti 
asymmetry A45Q while satisfying the constraints from the differential cross section. We further 
comment on the contribution to boosted top pair production. We consider several interesting 
limiting cases. The first is when there is no interference with the SM. Next the cases where the 
interference comes from only one operator are analyzed. Finally, we discuss the general scenario. 



5.1 No Interference: Cj^ = Cy = 

In models of heavy new physics which does not interfere with the SM, a simple relation between 
A45Q and A^7oo can be obtained from Eqs. (13) and (14): 

0.16cos2^A/'70o 



^450 



1 + 0.2AN; 



700 



(27) 
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Using the bound on A^7oo from Eq. (|4]) (which automatically satisfies the constraint from A^45o), we 
find A45Q < 0.07. From this we learn that the large ti asymmetry measured by CDF cannot be 
explained by a heavy new physics sector which does not interfere with the SM. This is consistent 
with the findings of |T9] , obtained directly from the data. Additionally, the maximal excess in the 
high-pr ti cross section is iVb ~ 2 . 



5.2 Vector Interference: c\ = and Cy ^ 



It is clear from Eq. ( 14 ) that Oy by itself does not contribute directly to A^gQ if = . However, 



taking Cy < can relax the constraints from A^45o and A^roo, thus allowing for a larger contribution 



to A^^Q from other operators, i.e. R ^ 0. Substituting R for Njqq and \cy \ for N^^q via Eq. (13) 



into the expression for A^^^q of Eq. ( 14 ) yields 



:i + Ar45o)xA 



COS 20 1.65 A^, 



700 



6.15 A". 



450 



19.67+ 14.65V1.8 + A", 



450 



0.24 A"- 



700 



(2^ 



Plugging in the constraints from Eq. (|4]), we find that the maximal value for A^^^q is 0.26 , which is 
1.3 (T below the mean value in Eq. ([2]). Regarding the boosted top cross section, an excess as large 
as A^b ~ 4 can be obtained. 
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Figure 1: The observables under consideration presented in the R — c\ plane for Cy = 6 = 0: The 
solid curve describes the mean value of A^^^q from Eq. ([2]), while the shaded region corresponds to 
the la range. The red-shaded region is the overlap of the latter with the la constraints on A450 
and A'700 in Eq. (Q. 



5.3 Axial Interference: c^j^O and Cy = 

Ref. [l2] showed that it is possible to explain the forward-backward asymmetry measurement with 
only 0\ . It is instructive to examine the addition of non-interfering operators. Fig. [l] shows the 
region in the parameter space of c\ and R satisfying Eqs. ^ and Interestingly, this region 
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is rather narrow, corresponding to ~ 2 and < 1 . Moreover, only the lower la range of A^jq 
can be accounted for in this case. (Note that the central value ^4450 ~ 0.4 requires a deviation of 
~ 1.8(7 in iVyoo , agreeing with [32] modulo the inclusion of l/A"^ effects.) As concerns the high-p-r 
tt cross section, a maximal excess of A*";, ~ 2 can be achieved (with R = 0) along with the la range 
for A« . 

5.4 The General Case 

We now explore the general parameter space accounting for the observables at hand. In Fig. [2] we 
show the allowed region in the Cy — c\ plane for various values of R . We learn the following: 

• As i? grows, the allowed region becomes smaller, and the maximal possible value is i? ~ 3.1 . 

• The allowed range for the vector octet operator is — 2 < Cy < . 

• The allowed range for the axial octet contribution is 0.3 < < 3.3 . 




Figure 2: The observables under consideration presented in the Cy — c\ plane: Each region corre- 
sponds to the overlap of the la ranges for A^^q , A450 and A'yoo in Eqs. ^ and Q, for different 
values of R . 



We conclude that: 

• In order to explain the measurement of the tt forward-backward asymmetry for Mu > 
450 GeV within la, a minimal contribution of the operator 0\ is necessary, ~ 0.3 . 
(This point in the parameter space corresponds to Cy ~ —1.9 and i? ~ 3.1 .) 

• The maximal enhancement of the boosted top pair cross section is A^t ~ 4 . Interestingly, 
this is consistent with within la . 

• Accounting for the high mass ti forward-backward asymmetry within la dictates a minimal 
excess of A^;, ~ 0.5 . 
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• Restricting the parameter space to include only an operator of definite chirality for each color 
structure, the maximal A^gQ is ~ 0.1 . An excess of Ai'f, ~ 4 can still be obtained. 



6 Predictions for Near-Future Measurements 



Thus far we focused on existing data from the Tevatron. This data could have interesting implica- 
tions on future Tevatron and LHC measurements. To illustrate this, we consider the tt differential 
cross section within our framework. We emphasize that our results below are general and include 
in particular the contributions from scalar and tensor operators. Fig. |3] depicts the Mu distribution 
at the Tevatron and LHC. The plotted regions correspond to the predicted enhancement relative 
to the SM, defined by 



tot 



da^^/dMu ' 

scanning over the entire parameter space obeying the la constraints of Eqs. ^ and (|4]). 



(29) 
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Figure 3: The ratio between the total and SM differential cross sections of top pair production 
as a function of Mti at the Tevatron (left) and the LHC at 7 TeV (right), calculated at leading 
order. The upper shaded regions correspond to the mean value of A^^q ~ 0.4 , scanning over the 
allowed range for A^45o and A^7oo defined in Eq. (|4]). The lower shaded regions correspond to the 
lower la range of A^I^q . The thick black curves at the bottom of the shaded regions correspond to 
i? = 0, 4 ~ -0.85, 4 ~ 1.7. 



We learn that if the latest high mass ti forward-backward asymmetry measurement persists 
and is accounted for by heavy new physics, then a significant enhancement of the ti differential 
cross section compared to the SM is expected at both the Tevatron and LHC. Specifically, at 
Mti ~ 1 TeV, a minimal factor of two enhancement is expected at the Tevatron. Similarly, at 
Mti ~ 1.5 TeV, the LHC should find a ti cross section of at least a factor of three higher than 
within the SM. In both cases the minimal enhancement is obtained for i? = 0, namely when only 
operators interfering with the SM are present. The combination that minimizes the enhancement, 
as described by the thick black curves at the bottom of the shaded regions of Fig. [3| does not 
depend on Mti and is given by Cy ~ —0.85, c\ ~ 1.7. To summarize: 

A^tot {Mu = 1 TeV) > 2 at the Tevatron , 
NtotiMtt = 1.5 TeV) > 3 at the LHC with = 7TeV . ^ ' 
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7 Outlook 



We have performed a model independent analysis regarding the tt forward-backward asymmetry, 
assuming heavy new physics. Any corresponding high scale new physics model can be mapped to 
our formalism to obtain constraints and predictions. We find a robust prediction in the form of 
enhancement in hard top physics at the Tevatron and LHC. The observation of such an enhance- 
ment would be exciting, and our analysis would assist in interpreting the signal and extracting 
microscopic information on the underlying physics. An equally intriguing possibility would be the 
absence of such an enhancement, assuming the asymmetry is established. Our findings would then 
imply the presence of sub-TeV new physics. Consequently, the new physics search strategy should 
be modified to include precision analysis in the 100-1000 GeV energy regime. 

Acknowledgments 

We thank Kfir Bhim, Alex Kagan, Zohar Komargodski, Seung Lee, Yosef Nir and Michele Papucci 
for useful discussions. GP is the Shlomo and Michla Tomarin career development chair and is sup- 
ported by the Israel Science Foundation (grant #1087/09), EU-FP7 Marie Curie, IRC fellowship, 
Minerva and G.I. P., the German- Israeh Foundations, and the Peter & Patricia Gruber Award. 

References 

[1] T. Aaltonen et al. [CDF Collaboration], Phys. Rev. D 83, 112003 (2011) [arXiv: 110 1.0034 
[hep-ex]]. 

[2] M. T. Bowen, S. D. Ellis and D. Rainwater, Phys. Rev. D 73, 014008 (2006) [arXiv:hep- 
ph/0509267]. 

[3] O. Antunano, J. H. Kuhn and G. Rodrigo, Phys. Rev. D 77, 014003 (2008) [arXiv:0709.1652 
[hep-ph]]. 

[4] L. G. Almeida, G. F. Sterman and W. Vogelsang, Phys. Rev. D 78, 014008 (2008) 
[arXiv:0805.1885 [hep-ph]]. 

[5] V. M. Abazov et al. [DO Collaboration], Phys. Rev. Lett. 100, 142002 (2008) [arXiv:0712.0851 
[hep-ex]]. 

[6] T. Aaltonen et al. [CDF CoUaboration] , Phys. Rev. Lett. 101, 202001 (2008) [arXiv: 0806. 2472 
[hep-ex]]. 

[7] CDF Collaboration, CDF note 10436, March 10, 2011. 

[8] CDF Collaboration, CDF note 10199, January 18, 2011; T. Aaltonen et al. [CDF Collabora- 
tion], arXiv:1106.5952 [hep-ex]. 

[9] CDF Collaboration, CDF note 10234, January 5, 2011. 

[10] K. Cheung and T. C. Yuan, Phys. Rev. D 83, 074006 (2011) [arXiv: 1101. 1445 [hep-ph]]. 

[11] C. Delaunay, O. Gedaha, S. J. Lee, G. Perez and E. Ponton, arXiv:1101.2902 [hep-ph]. 



10 



[12] J. Cao, L. Wang, L. Wu and J. M. Yang, arXiv:1101.4456 [hep-ph]. 

[13] Y. Bai, J. L. Hewett, J. Kaplan and T. G. Rizzo, JHEP 1103, 003 (2011) [arXiv: 1101.5203 
[hep-ph]]. 

[14] J. Shelton and K. M. Zurek, Phys. Rev. D 83, 091701 (2011) [arXiv:1101.5392 [hep-ph]]. 

[15] E. L. Berger, Q. H. Cao, C. R. Chen, C. S. Li and H. Zhang, Phys. Rev. Lett. 106, 201801 
(2011) [arXiv:1101.5625 [hep-ph]]. 

[16] M. L Gresham, L W. Kim and K. M. Zurek, arXiv:1102.0018 [hep-ph]. 

[17] V. Barger, W. Y. Keung and C. T. Yu, Phys. Lett. B 698, 243 (2011) [arXiv: 1102.0279 
[hep-ph]]. 

[18] B. Bhattacherjee, S. S. Biswal and D. Ghosh, Phys. Rev. D 83, 091501 (2011) [arXiv:1102.0545 
[hep-ph]]. 

[19] B. Grinstcin, A. L. Kagan, M. Trott and J. Zupan, Phys. Rev. Lett. 107, 012002 (2011) 
[arXiv: 1102.3374 [hep-ph]]. 

[20] K. M. Patel and P. Sharma, JHEP 1104, 085 (2011) [arXiv:1102.4736 [hep-ph]]. 

[21] G. Isidori and J. F. Kamenik, Phys. Lett. B 700, 145 (2011) [arXiv: 1103.0016 [hep-ph]]. 
[22] A. R. Zerwekh, arXiv:1103.0956 [hep-ph]. 

[23] E. R. Barreto, Y. A. Coutinho and J. Sa Borges, Phys. Rev. D 83, 054006 (2011) 
[arXiv:1103.1266 [hep-ph]]. 

[24] D. W. Jung, P. Ko, J. S. Lee and S. h. Nam, Phys. Lett. B 691, 238 (2010) [arXiv:0912.1105 
[hep-ph]]. 

[25] C. Degrande, J. M. Gerard, C. Grojean, F. Maltoni and G. Servant, JHEP 1103, 125 (2011) 
[arXiv:1010.6304 [hep-ph]]. 

[26] D. W. Jung, P. Ko and J. S. Lee, Phys. Lett. B 701, 248 (2011) [arXiv: 1011.5976 [hep-ph]]. 

[27] P. Ferrario and G. Rodrigo, Phys. Rev. D 78, 094018 (2008) [arXiv:0809.3354 [hep-ph]]. 

[28] P. Ferrario and G. Rodrigo, Phys. Rev. D 80, 051701 (2009) [arXiv:0906.5541 [hep-ph]]. 

[29] S. Jung, H. Murayama, A. Pierce and J. D. Wells, Phys. Rev. D 81, 015004 (2010) 
[arXiv:0907.4112 [hep-ph]]. 

[30] K. Cheung, W. Y. Keung and T. C. Yuan, Phys. Lett. B 682, 287 (2009) [arXiv:0908.2589 
[hep-ph]]. 

[31] P. H. Frampton, J. Shu and K. Wang, Phys. Lett. B 683, 294 (2010) [arXiv:0911.2955 [hep- 
ph]]. 

[32] J. Shu, T. M. P. Tait and K. Wang, Phys. Rev. D 81, 034012 (2010) [arXiv:0911.3237 [hep-ph]]. 



11 



[33] A. Arhrib, R. Benbrik and C. H. Chen, Phys. Rev. D 82, 034034 (2010) [arXiv:0911.4875 
[hep-ph]]. 

[34] I. Dorsner, S. Fajfer, J. F. Kamenik and N. Kosnik, Phys. Rev. D 81, 055009 (2010) 
[arXiv:0912.0972 [hep-ph]]. 

[35] J. Cao, Z. Heng, L. Wu and J. M. Yang, Phys. Rev. D 81, 014016 (2010) [arXiv:0912.1447 
[hep-ph]]. 

[36] V. Barger, W. Y. Keung and C. T. Yu, Phys. Rev. D 81, 113009 (2010) [arXiv: 1002. 1048 
[hep-ph]]. 

[37] Q. H. Cao, D. McKeen, J. L. Rosner, G. Shaughnessy and C. E. M. Wagner, Phys. Rev. D 
81, 114004 (2010) [arXiv:1003.3461 [hcp-ph]]. 

[38] B. Xiao, Y. k. Wang and S. h. Zhu, Phys. Rev. D 82, 034026 (2010) [arXiv: 1006.2510 [hep-ph]]. 

[39] R. S. Chivukula, E. H. Simmons and C. P. Yuan, Phys. Rev. D 82, 094009 (2010) 
[arXiv:1007.0260 [hep-ph]]. 

[40] C. H. Chen, G. Cvetic and C. S. Kim, Phys. Lett. B 694, 393 (2011) [arXiv:1009.4165 [hep- 
ph]]. 

[41] E. Alvarez, L. Da Rold and A. Szynkman, JHEP 1105, 070 (2011) [arXiv:1011.6557 [hep-ph]]. 

[42] K. Blum et ai, Phys. Lett. B 702, 364 (2011) [arXiv: 1102.3133 [hep-ph]]. 

[43] T. Aaltonen et al. [CDF Collaboration], Phys. Rev. Lett. 102, 222003 (2009) [arXiv: 0903. 2850 
[hep-ex]]. 

[44] V. Ahrens, A. Ferroglia, M. Neubert, B. D. Pecjak and L. L. Yang, JHEP 1009, 097 (2010) 
[arXiv: 1003.5827 [hep-ph]]; arXiv: 1103.0550 [hcp-ph]. 

[45] M. Cacciari, S. Frixione, M. L. Mangano, P. Nason and G. Ridolfi, JHEP 0809, 127 (2008) 
[arXiv:0804.2800 [hep-ph]]. 

[46] N. Kidonakis, Phys. Rev. D 82, 114030 (2010) [arXiv: 1009.4935 [hep-ph]]. 

[47] U. Langenfeld, S. Moch and P. Uwer, Phys. Rev. D 80, 054009 (2009) [arXiv:0906.5273 [hep- 
ph]]. 

[48] Y. Eshel, O. Gedaha, G. Perez and Y. Soreq, Phys. Rev. D 84, 011505(R) (2011) 
[arXiv: 1101. 2898 [hep-ph]]. 

[49] N. Kidonakis and R. Vogt, Phys. Rev. D 68, 114014 (2003) [arXiv:hep-ph/0308222]. 

[50] V. Khachatryan et al. [CMS Collaboration], Phys. Rev. Lett. 106, 201804 (2011). 
[arXiv: 1102.2020 [hep-ex]]. 

[51] G. Aad et al. [ATLAS Collaboration], New J. Phys. 13, 053044 (2011). [arXiv:1103.3864 

[hep-cx]]. 

[52] A. D. Martin, W. J. Stirling, R. S. Thorne and G. Watt, Eur. Phys. J. C 63, 189 (2009) 
[arXiv:0901.0002 [hep-ph]]. 



12 



